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Abstract 

Experiments searching for rare events, such as neutrinoless double beta decay and interac- 
tions of dark matter candidates, require extremely low levels of background. When these ex- 
periments are performed using macro-bolometers, radioactive contamination near the surfaces is 
of particular concern. For a bolometric neutrinoless double beta decay experiment, it can cause 
counts in the spectral region where the signal is expected, while for a dark matter experiment 
which exploits ionization signals for particle identification, it originates an incomplete charge 
collection simulating a nuclear recoil. In order to control the effects of surface contamination, 
we developed a novel technique that uses composite macro-bolometers to identify energy depo- 
sitions that occur close to the surfaces of materials immediately surrounding the detector. The 
composite macro-bolometer proposed and studied here consists of a main energy absorber that is 
thermally coupled to and entirely surrounded by thin absorbers that act as active shields. Surface 
energy depositions can be rejected by the analysis of simultaneous signals in the main absorber 
and the shields. In this paper, we describe a full thermal model and experimental results for three 
prototype detectors. The detectors consist of Ge, Si, or Te02 thin absorbers as active shields, all 
with Te02 crystals as main absorbers. In all cases, the surface event rejection capability is clearly 
demonstrated. In addition, simulations and preliminary results show that it is possible to detect 
energy depositions that occurred on the shields without separate readout channels for them. The 
energy depositions in the shields are distinguished from those in the main absorber through pulse 
shape discrimination. This simplification makes this technique a viable method for the rejection 
of surface energy depositions in next-generation bolometric double beta decay searches, such as 
possible extensions or upgrades of the CUORE experiment. 
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1. Introduction and motivation 

The searches for the neutrino mass and dark matter are at present two of the most relevant 
and exciting fields in cosmology and particle physics. Experiments that search for neutrinoless 
double beta decay (0y/3j3) [1] or the nuclear recoils induced by a WIMP [2] require the detection 
of very rare events. Although the signals are expected in two very diff'erent energy regions, 
~ MeV for the former and < 100 keV for the latter, the experimental approaches to these searches 
often share common technological challenges. In both cases, bolometers are used as detectors 
for many sensitive experiments, such as EDELWEISS [3], CRESST [4] and CDMS [5] for dark 
matter searches or Cuoricino [6] and CUORE [7] for neutrinoless double beta decay. 

Bolometers are phonon-mediated particle detectors operated at low temperatures [12]. These 
devices are capable of obtaining both higher energy resolutions and lower energy thresholds than 
conventional detectors. In addition, they can be fabricated from a wide variety of materials, 
allowing flexibility for experiments that require the detectors to contain particular atomic or 
nuclear species. If other excitations (such as ionization charge carriers or scintillation phonons) 
are exploited in addition to phonons, bolometers can discriminate nuclear recoils from electron 
recoils, or a particles from particles and y rays. 

Bolometer-based OvjSlS and dark matter experiments require extremely low levels of radioac- 
tive background, since they search for very rare events. For difl'erent reasons in the two cases, 
surface radioactive contamination is of particular concern. In OyjSfi searches, a's and jS's of su- 
perficial origin can lose part of their energy in a few microns and generate a continuum in the 
spectral region where the Ov^jS signal is expected. In dark matter searches, electron recoil at the 
detector surface (from any type of ionizing interactions) can result in incomplete charge collec- 
tion, mimicking a nuclear recoil. Therefore, a bolometer capable of tagging surface events would 
be a powerful tool to improve the signal-to-background ratio in both types of experiments. 

In monolithic bolometers it is possible to obtain some spatial information by detecting phonons 
before they thermalize. This is used in the dark matter experiments CDMS and EDELWEISS [13, 
14]. Recently, it was shown that this technique can be extended to bolometric detectors for OyjSjS 
as well, though incomplete thermalization leads to degradation in energy resolution [15, 16]. For 
dobule beta decay experiments, it is important to have good energy resolution to improve the 
signal-to-background ratio. In this case, it is preferable to operate the bolometers in calorimetric 
mode, where the phonon thermalization from an event is almost complete and the deposited en- 
ergy is fully converted to heat and measured as a temperature rise. However, this mode loses all 
spatial information about the event, which makes rejection of surface background events more 
difficult. 

The aim of this study is to demonstrate that the lack of spatial resolution intrinsic to low tem- 
perature calorimeters may be partially overcome by utilizing a composite bolometer consisting 
of a main energy absorber that is thermally coupled to and surrounded by thin absorbers working 
as auxiliary bolometers and active shields (Fig. 1). We call this composite bolometer a "surface 
sensitive bolometer" or SSB. 

The motivation of this work originates in physics experiments being carried out at the under- 
ground Gran Sasso National Laboratory (LNGS) in Italy searching for neutrinoless double beta 
decay of ^^^Te and named Cuoricino [6] and CUORE [7]. Cuoricino ended in June 2008 and 
was constituted by 62 Te02 crystals of two sizes, ~ 330 g and ~ 790 g, with a total of 40.7 kg of 
Te02. CUORE is a next generation experiment currently under construction, and will have 988 
crystals for a total of 741 kg of Te02. Cuoricino was operated at 10 mK, and CUORE is designed 
to operate at similar temperatures. These experiments look for a peak in the energy spectrum at 
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Figure 1 : Schematic 3D view of a surface sensitive bolometer. The main absorber is shown in dark gray (blue) while the 
thin auxiliary bolometers are in light gray (red) and are thermally connected to the main absorber. The small elements 
on each absorber surface are the temperature sensors. 

2527 keV [10, 11]. The observation of Qvpp peak would imply that neutrinos are self-conjugate 
particles and enables a sensitive measurement of their mass scale. These experiments are looking 
for events of less than 10"^ counts/y/mole and it is critical to control the background generated 
by residual radioactive impurities both inside and near the detector. A background of at least 0.01 
counts/keV/kg/y near the energy region of interest is required for CUORE to reach its desired 
sensitivity. 

Simulations based on the background levels observed in Cuoricino and measurements of 
radioactive contaminants with germanium detector gamma counting and other techniques show 
that ~ 60% of the contribution to the radioactive background at the energy of interest can be 
identified as partial energy deposition of a and yS-particles emitted from either the surface of 
the detectors or the materials that surround them [6]. These particles release only a fraction 
of their energy into the bolometer, and produce a roughly continuous spectrum which extends 
from the full energy of the decay to zero. With respect to Cuoricino, a substantial improvement 
was obtained in R&D tests performed in view of CUORE [8]. In particular, the contribution 
coming from Te02 crystal surfaces was reduced thanks to a new polishing procedure [9] and is 
now compatible with the aforementioned target of 0.01 counts/keV/kg/y. On the contrary, the 
background coming from the inert material surfaces around the detector, in particular the Cu 
holding structure, has not been reduced yet at the desired level. The work described here aims at 
decreasing this background by developing detectors that are capable of identifying these surface 
events [17, 18]. 

The present paper is organized as follows. In Section 2 the basic principles and structure 
of the proposed detectors are described. Afterwards (Section 3), we present a thermal model of 
the SSB that we use to simulate its performance. The rest of the paper describes our prototype 
detectors equipped with SSBs and their experimental characterization. After a description of the 
setup in Section 4, the experimental data are shown and discussed in Section 5 with a presentation 
of two approaches for the identification of surface events and a preliminary quantification of the 
rejection power. Section 6 is dedicated to a further approach that promises to tag surface events 
with a passive use of the shields. Preliminary encouraging results and simulations are presented. 
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Conclusions and prospects are discussed in Section 7. 

2. Basic structure and operation principle of Surface Sensitive Bolometers 

A bolometer is a solid-state particle detector composed of two elements: an energy absorber, 
where the energy of the particles is deposited and converted into phonons, and a phonon sensor 
attached to the absorber. The particle energy is initially stored in the form of optical and Debye- 
energy phonons. If the absorber is kept at very low temperatures (below 100 mK), the phonons 
created by the particle interaction have much higher energies than the phonon thermal bath and, 
in bolometers, represent the elementary excess excitations which mediate the particle detection. 
For macro-bolometers with dimensions of a few centimeters as those considered here, the phonon 
thermalization time is typically of the order of few microseconds. If either the phonon response or 
the phonon transmission to the sensor is slower, then the deposited energy has enough time to be 
converted to heat and the phonon sensor interprets a particle energy deposition as a temperature 
rise, i.e. the phonon sensor effectively acts as a thermometer. In the present work, phonon 
sensors are glued with thick (up to 50 jim) epoxy layers at the energy absorber, introducing an 
intrinsically slow (~ 1-10 ms scale) transmission interface. The bolometer is then operated in a 
calorimetric mode which provides the optimal energy resolution. 

Resistive elements with a strong temperature dependence are often used as thermometers. 
There are two types of thermometers typically used for these experiments: semiconductor ther- 
mistors, like the ones used in this study, or transition edge sensors (TES), which are super- 
conducting films kept at their critical temperature. To minimize heat capacity, bolometers are 
typically operated at cryogenic temperatures below 100 mK (in some cases, less than 15 mK). 
Dielectric diamagnetic crystals have low heat capacities and are commonly used as energy ab- 
sorbers. With such devices, energy resolutions as low as 5 eV have been achieved for X-rays [19] 
and 5 keV for a-particles [20] ; y-rays can be detected with resolutions comparable to those ob- 
tained with the best germanium diodes [21]. 

In calorimetric operation, no diff'erence is expected for the detector response to interactions 
of diff'erent ionizing particles, in particular electrons and Of-particles. However, if a relevant frac- 
tion of athermal phonons contributes to the signal formation, in principle diff'erences between 
electrons and a-particles can arise, since the diff'erent track structure and ionization density could 
influence the heating process. This aspect was studied experimentally by comparing the response 
of Te02 bolometers to y-rays and a-particles [22] . It was found that the amplitudes of the signals 
provided by a glued semiconductor thermistor are the same within 2% for the same energy de- 
position of y-rays and a-particles. No appreciable difl'erence in the signal shapes was observed. 
On the basis of this result, the detector responses to electrons and a-particles will be considered 
to be the same in the following discussion. 

The basic idea of the SSB is to surround the main absorber with thin auxiliary absorbers 
which function as active shields and are operated as bolometers. In this study the main absorber 
and the shields are thermally connected, and a near 4;r coverage from external charged particles 
can be achieved. The SSB is designed to be able to distinguish among three types of events: 
those induced by particles originating outside of the entire detector and that are stopped in the 
shields depositing only in them their residual energy, those induced by particles originating on the 
surface of the main or auxiliary absorber with energy deposited in both, and those with energy 
deposited in the bulk of the main absorber. We refer to the first case as "surface events", the 
second one as "mixed events" and the third one as "bulk events". 
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The main absorber and the active shields may be made of the same or different materials 
depending on purity requirements and appropriate thermal properties. Because the main interest 
of this study is to reduce the background for CUORE, the main absorber proposed here and 
studied in the model is made of Te02 and is cubic or rectangular in shape. Each of its six surfaces 
will be covered with an active shield, the area of which is similar to that of the corresponding face. 
In order to not significantly alter the detector geometry and space occupancy, the shields have a 
thickness of a fraction of a millimeter, about two orders of magnitude less than the main absorber 
dimension. Active shields of Si, Ge, and Te02, are considered. The choice of these materials, 
as discussed more extensively in Section 4.1, was mainly driven by reasons of radiopurity (Ge), 
availability and cost (Si), and thermal compatibility (Te02).In the current design, the shields are 
thermally and mechanically coupled to the main absorber by epoxy beads of thickness less than 
~ 50 Jim. Other coupling methods between shields and main absorber can be devised, aiming at 
a better control of the thermal conductance between these elements. One of them is implemented 
in this work and described in Section 4.2. 

The thermal pulses from energy deposition in the absorbers are read out by thermistors ther- 
mally coupled to the absorbers by epoxy beads. The bolometers described in this work use 
semiconductor thermistors operated in the variable range hopping (VRH) conduction mode [23]. 
These thermistors consist of neutron transmutation doped (NTD) Ge single crystals, with a mass 
of the order of ~ 10 mg. The doping technique [24] allows to achieve uniform doping in the 
whole thermistor volume and to accurately tune the net dopant concentration, corresponding to 
the desired temperature dependance of the resistivity at low temperatures. 

The main absorber and each of the six shields have their own thermistor, and each absorber- 
thermistor pair can be regarded as an individual detector. The bias on the thermistors results in 
a static heat flow and therefore a temperature diff'erence between the thermistor and the absorber 
appears. These difl'erences are generally much smaller than the absolute temperatures of the 
detector components. The main absorbers are thermally linked to a copper frame acting as a 
heat sink through a set of PTFE (Polytetrafluoroethylene) stand-off's, and the copper frame is 
thermally coupled to the mixing chamber, the coldest point of a dilution refrigerator [25]. 

The main absorber, epoxy beads, and thermistors are the same as those used in Cuoricino, 
and their thermal properties measured for Cuoricino have been used as input for the simulations 
described in Section 3. 

Charged particles from materials outside the SSB, such as a-particles, will be stopped and 
tagged by one of the active shields (surface events). They release part to most of their energy in 
the shield, but raise the temperature of all the detector elements as they are, in fact, all thermally 
connected. Because of the small heat capacity of the shield due to its small mass, the signal 
read by its thermistor will have a higher amplitude and faster rise time than the signal read by 
the thermistor attached to the main absorber. If, on the other hand, an energy deposition occurs 
inside the main absorber (bulk event), all of the thermistors will read pulses with comparable 
amplitudes and rise times. The rise time of the temperature pulse in the active shield thermistor 
is much slower when energy is released in the bulk of the main absorber as opposed to when it 
is released in the shield. It is therefore possible to separate bulk and surface events by compar- 
ing amplitude and shape of pulses among the different thermistors. No degradation in energy 
resolution of the main absorber is expected since the device is still operating in the mode where 
phonon thermalization is nearly complete. 

A specific discussion is required for a-particles that originate from surface contamination of 
the main absorber or of the shield itself. As already pointed out, one of the largest backgrounds 
for neutrinoless double beta decay experiments like Cuoricino comes from surface contaminants 
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that emit a-particles. The a-particles can mimick a OvySyS event if they deposit a part of their 
energy elsewhere and about 2530 keV in the main absorber, where the double beta decay peak is 
expected. By using the active shields discussed here, a's originating from outside the detectors 
or from the external surfaces of the shield can be distinguished from the bulk events in the main 
absorbers. As for a's escaping from the internal surfaces of the shield or from the energy absorber 
surface, and delivering in the main crystal about 2.5 MeV (close to the double beta decay signal), 
they will surely deliver at least 1.5 MeV in the shield (their remaining energy), since natural a's 
have energies higher than ~4 MeV. These mixed events can be easily recognized by the large 
fast pulses associated that are characteristic of the energy deposited in the shields. We point out 
however that, as discussed in Section 1 , the background originating from the Te02 surface is now 
under control thanks to appropriate cleaning methods. 



3. Thermal model and expected performances of Surface Sensitive Bolometers 

For a thermal detector without the active shields, the thermal network consists of three ther- 
mal nodes: the main energy absorber (here a Te02 crystal), the thermistor lattice, and the ther- 
mistor electrons (Fig. 2a). The thermistor is split into two elements because it has been observed 
that in these devices the lattice phonon system and the conduction electron system behave as 
separate thermal stages that reach different temperatures and are connected by a finite thermal 
conductance, which is internal to the thermistor itself. Such decoupling is due to a non-ohmic 
eff'ect explained by the so-called "hot electron model" [26]. The electron-phonon thermal con- 
ductance of thermistors of the same type as those used in this work was measured. The methods 
and the results of these measurements are reported elsewhere [27, 28]. It has also been shown 
experimentally that there is no direct connection between the thermistor electron system and the 
heat sink [27]. 

In Fig. 2a, the presence of heat conductors in the system is shown. As described in section 2, 
the heat conductors between thermistors and absorbers are epoxy beads, while PTFE mechanical 
supports for the main absorber act as heat conductors between the main absorber and the heat 
sink. The thermistor read-out wires connect the thermistor lattice system to the heat sink. 

The thermal model for the SSB is an extension of the three-node model and requires three 
additional nodes for each active shield used in the system [29]. For simplicity, the network 
adopted here includes the absorber and thermistor for the main detector and one active shield 
(Fig. 2b). In general, the heat conductor between the main absorber and the active shield is again 
a set of epoxy beads, although other solutions can be devised (see Section 4.2). 

In a static condition with external power tt/ injected into node / and power Ptj flowing between 
nodes / and j, the power balance is given by: 
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''i-Tj^^j = ^ /=1,2,...,6 (1) 

where ^ 

Ptj= f Gtj(T)dT (2) 

Here, Ti is the equilibrium temperature of node /, and Gij is the temperature-dependent thermal 
conductance between nodes / and j, with Gij = Gji and Gu = 0. Gij = whenever there is 
no physical direct connection between nodes / and j. The node 7 = corresponds to the heat 
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Figure 2: (a) Three-node thermal model for a bolometer, consisting of the main absorber node and the two thermistor 
nodes. The thermistor is separated into electron and phonon systems, connected by the electron-phonon thermal con- 
ductance, (b) Thermal network used for the case of a single-shield SSB. The presence of the shield adds three nodes, 
corresponding to those of the shield-thermistor system. In our experiment, epoxy beads are normally used to obtain 
heat conduction between the main absorber and the shield (Gmain-shieid) and between the thermistors and the absorbers 
(Gabs-ph, ^abs-p})' addition, PTFE holders link the main absorber to the heat bath {G abs-b) and gold wires provide 
electric read-out and thermal connection between thermistor lattice and heat bath (Gph-b, ^ph-b^' ^ ^^^^ 
capacities of each element of the main and of the shield bolometer respectively. 



bath. 7ii also includes the parasitic power dissipated in the absorber (mainly due to mechanical 
frictions between the main absorber and the PTFE supports) and in the electron system (due 
probably to parasitic currents induced by EM interferences). The thermistors are biased with 
constant current, and the detector as a whole reaches thermal equilibrium with the six thermal 
stages at different temperatures. The electric power dissipated in the thermistors from the bias 
current is also included in the tt/ for the electron system of the thermistor. For well fabricated 
and shielded bolometers, the dominant power dissipation is ascribable to the bias current. 

The dynamic behavior of the detector must be determined to evaluate the detector response 
to energy deposited by particles interacting either in the main absorber or in the active shield. 
The corresponding differential equation for node / is 

6 

Ci(Ti) . ti = ni-Y, Pij / = 1 , 2, . . . , 6 (3) 

where C/ is the heat capacity of node /. An explicit dependence of the term tt/ on the node tem- 
perature Ti{t) has to be introduced for thermistor nodes connected to a bias circuit to account for 
a characteristic mechanism known as electrothermal feedback. This additional effect is caused 
by power dissipation by the bias current which raises the thermistor temperature and acts back 
on its resistance until an equilibrium is reached. This mechanism, responsible for the deviation 
of thermistor voltage-current curve from linearity, is included in the calculations. 

The time evolution of a heat pulse is computed by solving a set of differential equations with 
the initial conditions provided by the solution of the static problem and by the temperature rise 
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Table 1: Detector parameters C, G for Eqs.(l-3) used in the simulations. The simulated SSB, represented by the thermal 
network in Fig. 2b, has a Te02 cubic main absorber with 5 cm side length and a single silicon active shield with 5x5 cm^ 
surface and 300 jim thickness. T is given in Kelvin. The heat sink is assumed to be at the base temperature T}j=9 mK. 
The Rq, To parameters are those used in Eq.(4). 
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due to the energy deposited by the absorbed particle. Temperature pulses on the thermistor- 
electron systems of both main absorber and active shield are then converted into voltage pulses 
using the resistance vs. temperature behavior of the thermistors parameterized by 

R(T) = Ro exp[(To/Tf'] (4) 

as predicted by VRH with Coulomb gap [23]. The conversion is possible because of the constant 
current / flowing in the thermistor. In fact, a resistance variation AR due to a temperature change 
Ar of the electron system results in a voltage variation AV=I • AR=I • (dR/dT) • AT. In the 
experiements, typical pulse amplitudes are of the order of several tens of microvolts per energy 
deposition of 1 MeV. These voltage pulses are amplified to a range of a few volts by a low noise 
voltage amplifier. 

Numerical codes were developed for calculations of both static and dynamic behaviors. Val- 
ues for thermal conductances and heat capacities were obtained experimentally from previous 
measurements with detectors similar to those used in Cuoricino [27] . A set of simulations was 
performed to determine the dependance of the discrimination efficiency on several experimental 
parameters. The simulations are based on the same detector scheme: one Te02 main absorber, 
and one shielding element. A pair of thermistors with the same Rq and Tq were used for the main 
absorber and its shield. 

We have evaluated the behavior of a detector with the configuration adopted for Cuoricino 
and give the relevant thermal properties in Table 1 . We have simulated a Si active shield attached 
to the main Te02 crystal with a surface area of 5x5 cm^ and a thickness of 300 yum. The Cuoricino 
Te02 crystal is a cube with 5 cm side length. The same size is foreseen for the CUORE detectors. 

Fig. 3 shows three pairs of simulated pulses. They are read by thermistors on both the main 
absorber and the active shield simultaneously, following a specified energy deposition in the 
detector. Fig. 3(a) shows simultaneous pulses for a 2.5 MeV energy deposition in the Te02 crystal 
(bulk event). Fig. 3(b) shows pulses from an energy deposition of 2.5 MeV shared in equal parts 
between shield and absorber (mixed event), and Fig. 3(c) shows pulses from a 2.5 MeV energy 
deposition in the shield only (surface event). The first pulse corresponds to the signal, while 
the second and third ones to background. The simulated pulses from the main absorber are 
consistent with those observed in Cuoricino both in amplitude and time structure: the amplitude 
is ~ 140 yuV/MeV, the rise time (10%-90%) is ~ 50 ms, and the decay time (90%-30%) is 
-500 ms [29]. 
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Figure 3: (a) Pair of simulated pulses generated by an energy deposition of 2.5 MeV in the main absorber of the SSB. 
The pulses, which are read by the thermistor of the main absorber and by the thermistor of the active shield, present 
comparable amplitudes and shapes. The temperature increases corresponding to the voltage pulses are 5 • 10-2 mK for 
the main absorber and 3 • 10"^ mK for the shield, (b) Pair of simulated pulses generated by a simultaneous energy 
deposition of 1.25 MeV in the shield of the SSB and of 1.25 MeV in the main absorber, (c) Pair of simulated pulses 
generated by an energy deposition of 2.5 MeV in the shield of the SSB. Compared with the main absorber, the signal seen 
by the shield thermistor is higher and faster due to the shield lower heat capacity. The temperature rises corresponding 
to the voltage pulses are 3 • 10"^ mK for the main absorber and 1.9 mK for the shield. 

When the energy is deposited in the main absorber, the pulse decay time is ~ Cabs/Gabs-b, 
and is of the order of ~ 500 ms at 10 mK. (Here and in the following the symbols for the 
thermal parameters are the same used in Fig. 2.) The corresponding rise time is given by 
C parallel I Gph-eu whcrc C parallel = Cabs ' Cell {Cabs + Q/). This time constant Can be approximated 
by Cei/Gph-ei ~ 80 ms. If the energy is released in the active shield, Cabs is replaced by the heat 
capacity of the shield C\bs, two orders of magnitude lower than that of the main absorber. The 
ratio of the main-absorber to shield is {CabslC\bs) ' (C\bs + Cei)/(Cabs + Cei), which can be 
approximated, considering the values of the heat capacities in our model, with CeilC\bs'- this is 
of the order of 10^ at 10 - 15 mK. These time constants, determined with simple approximations, 
are in reasonable agreement with the detailed simulation results shown in Fig. 3. 

Information on the site of an energy deposition is given by comparing the voltage pulse am- 
plitude read by the thermistor on the main absorber with the simultaneous voltage pulse ampli- 
tude read by the thermistor on the active shield. One can then draw a scatter plot of shield pulse 
amplitudes on the y-axis versus simultaneous main pulse amplitudes on the x-axis. The relation 
between the amplitudes is shown by the curves in Fig. 4. The simulation shows different classes 
of pulses: the upper curve results from energy depositions inside the shield only, the lower one 
from energy depositions inside the main absorber. The intermediate curve represents events due 
to simultaneous energy depositions of equal energy both in the shield and in the main absorber 
(mixed events). The shielding element can reach saturation in the thermistor due to the low heat 
capacity of the shield. This simulation shows that in principle the origin of energy depositions 
can be determined through a scatter plot analysis. 

Simulations also show the possibility of discriminating between shield and main-absorber 
events using the information from the shield thermistor alone, without the need for a separate 
thermistor on the main absorber. Energy depositions in the shield give rise to a higher and faster 
signal than those happening in the bulk. 
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Amplitude V^on main absorber thermistor [mV] 

Figure 4: Relation between peak values of the voltage pulses from main absorber and from the shield. The upper curve 
results from energy depositions in the shield only, the lower one from energy depositions in the main absorber. The 
intermediate curve represents events due to simultaneous energy depositions of equal energy both in the shield and in the 
main absorber. The energy range of this simulation is from 100 keV to 9 MeV. As an example, points corresponding to 
an overall energy deposition of 2 MeV are circled. In all cases, a high-energy saturation of the curves is apparent, due 
to the non-linearity induced by the large temperature signals Ar. This effect is not expected in the experimental results 
here described, since in this case the much higher operation temperature (~25 mK versus ~ 10 mK) leads to much smaller 

Ar/r. 

4. Design and realization of the experimental setup 

Three prototype SSBs were assembled and tested. All prototypes used Te02 crystals as 
the main absorber. Three different materials, Ge, Si and Te02, were used as shield absorbers. 
The scope of this multiple choice is to show experimetally that the rejection of surface events 
can be efficiently achieved irrespectively of the shield composition, providing more degrees of 
freedom in the design of future devices. The three chosen materials have pros and cons, below 
summarized. 

4.1. The choice of the material for the shields 

The advantage of Ge and Si is related to the very well developed crystallization technology. 
Large mass crystals can be grown, and the achievement of thin wafers from them is a routine 
operation in semiconductor industry. In addition, sophisticated methods for the treatment of the 
surfaces (grinding, polishing, etching), particularly relevant in this context, are available and 
spreadly used. From the bolometric point of view, Ge is a commonly used absorber with high 
performance [3]. Si is also successfully employed in bolometric experiments [5], and, at least 
on the paper, it should provide higher and faster pulses than Ge in equal conditions, due to 
its higher Debye temperature (645 K versus 360 K). In practice, both materials are excellent 
for the operation of low temperature macro-calorimeters, and Ge is better suited to the use of 
NTD Ge sensors since the mechanical coupling between equal-material elements does not pose 
particular problems. From the point of view of intrinsic purity, Ge is by far the purest material 
achievable, thanks to the zone-refining crystallization method. This reflects also on the level of 
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the attainable radiopurity, which is excellent and well below 10"^^ g/g in U and Th concentration. 
Furthermore, the concentration of radioisotopes is well known thanks to the high sensitive Ge- 
based Ovpp experiments performed so far or currently in operation (IGEX, Heidelberg-Moscow 
and GERDA [1]). The intrinsic radioactivity of Si is not known with the same accuracy, however 
nothing prevents in principle from obtaining similar low levels. In addition, high purity Si can be 
produced at prices significantly lower. This element is relevant when planning experiments with 
hundreds of channels, as in the present case. 

However, a significant disadvantage of these two materials in the application here discussed 
has to be underlined. Both Ge and Si present isotropic thermal contraction due to their crystalline 
structure. On the contrary, the Te02 main absorber has equal thermal expansion coefficient along 
directions orthogonal to the [001] growth axis, but a diff'erent (and furthermore negative) one 
along it. This means that the mechanical coupling of a thin Ge or Si element, which is anyway 
difficult at surfaces orthogonal to the [001] axis because of the diff'erent expansion coefficients, 
becomes even problematic at those surfaces which are parallel to the growth axis. We experi- 
mented this problem in the present work, often observing a spontaneous detaching of the Ge or 
Si shields after thermal cycling. If these materials were chosen for the final application, a specific 
R&D activity has to be considered to design safely the mechanical and thermal coupling of the 
shields. 

The advantage of the Te02 shields is of course immediately related to the main (and only) 
disadvantage of the Ge and Si ones: they match exactly the thermal contraction features of 
the main absorber, after taking care of producing two types of shields (respectively orthogonal 
and parallel to the growth axis) and of coupling them to the appropriately corresponding main- 
absorber surfaces. From the point of view of the bulk radiopurity, we can consider crystalline 
Te02 as a safe and well known material (although not at the level of Ge and Si), thanks to the 
successfull Qvpp experiments performed up to now [1,6] and to the extensive study conducted 
by the CUORE collaboration on this subject [9]. The adequate level of the bolometric perfor- 
mance of Te02 is of course out of discussion. The issue here is related to the production of thin 
large-surface slabs with an appropriate sawing method, capable to introduce negligible surface 
radioactivity. Discussions with the Te02 crystal company selected by the CUORE collaboration^ 
and a preliminary production show that this operation is indeed possible. The fragility of the thin 
Te02 shields could also be a problem for an easy mechanical assembly of the detectors, but the 
experience gathered in this work is reassuring under this aspect. 

4.2. Description of the detectors 

Fig. 5 shows schematics of the three detectors and Table 2 shows their operating properties. 
Fig. 6 shows a photograph of the Si SSB. Each main absorber was sandwiched between two 
active shields, and the main absorber is thermally and mechanically coupled to the heat bath by 
four PTFE supports. The tests described here were conducted in the Cryogenic Laboratory at 
Universita dell'Insubria at Como, Italy. 

In order to avoid detector saturation from cosmic rays, 2 x 2 x 0.5 cm^ and 2 x 2 x 0.8 cm^ 
main absorbers (see Table 2) were used instead of the 5 x 5 x 5 cm^ cubes used in Cuoricino. 
Cuoricino is located underground at Laboratori Nazionali del Gran Sasso (LNGS), Italy, and 
is shielded from cosmic rays under 3500 m.w.e. of rocks, while the laboratory at Insubria is 
aboveground. 



^SICCAS, Shanghai Institute of Ceramics - Chinese Academy of Sciences - Shanghai, China 
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Figure 5: Schematic of the three SSB prototypes. Each detector consists of one main absorber and two active shields. 
The active shields are made of Ge, Si and Te02 respectively. 
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Figure 6: Photograph of a surface sensitive bolometer made of Te02 and Si active shields. The light grey rectangular box 
is the main absorber (20x20x 5 mm), held by four PTFE supports 5 mm wide. One of the two Si shields is visible (the 
other one is on the opposite face of the Te02 crystal). The Si shield is glued at the 20x20 mm face of the main absorber 
and appears as a square reflective thin layer with a 15x15 mm surface area. 
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Table 2: Main properties of the three surface sensitive bolometers used for the experimental tests: active shield material, 
volume and mass of the main absorber, volume, mass and heat capacity (computed at 25 mK) of the shield, thermistor 
volume, and resistance-temperature VRH parameters of the thermistor - see Eq. (4). 
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The shields are 0.5 mm thick for Ge and Te02 and 0.3 mm for Si, and thermally coupled to the 
Te02 main absorber by four epoxy beads of 1 mm diameter and 50 jim thickness. The mechanical 
mounting is designed to minimize the risk of the shields detaching due to the difference in thermal 
contractions at cryogenic temperatures. This is especially true for the Ge and Si, as discussed 
above. An assembly method involving relatively thick epoxy beads was chosen in the Ge and Si 
case because the glue disks proved to be sufficiently soft to accomodate the differential material 
contractions. The shields were glued at surfaces orthogonal to the [001] axis, for which the stress 
induced by the thermal contractions is less critical. 

For Te02 shields, the [001] axis of the main and shield absorbers must be aligned to match the 
thermal contraction when cooling to low temperatures, as discussed in the previous subsection. In 
this case, where no difference is expected between the thermal expansion coefficients of shields 
and main absorber, we adopted a different gluing method, which implies a more rigid adhesion. 
Four Ge stand-offs (1 mm^ in area and 50 jim thick) were inserted between the Te02 shields and 
the main absorber, with thermal coupling established by the same epoxy used in the Si and Ge 
case, but this time in the form of a thin veil. These stand-offs were introduced to provide a more 
reproducible thermal connection between the shields and main absorber, since the contact area 
is precisely defined by the stand-off geometry. As for Si and Ge shields, the amount of glue was 
carefully controlled. It was in this case much less, due the largely inferior thickness of the epoxy 
layer. We expect that the thermal conductance through the stand-offs, never measured, is of the 
same order of magnitude as that of the beads, since the interface glue-crystal should provide the 
dominant contribution. 

The main absorber and both shields are thermally coupled to their own NTD Ge thermistor. 
The thermal coupling between absorbers and thermistors is made by six epoxy beads (0.5 mm 
diameter and 50 jim thickness) for the main absorbers and Ge and Si shields and a single larger 
bead for the Te02 shields. When six beads are used, they are placed on the 3 x 1.5 mm ^ surface 
of the thermistor (see Table 2 for thermistor size). 

An external source of a-particles was used to test the SSB. The a source was a piece of copper 
strip implanted with ^^^Ra. ^^^Ra emits Of-particles with a half life of 3.66 days in equilibrium 
with its a and jS emitting daughters. The main a lines are at 5.68, 6.29, 6.78 and 8.78 MeV. 
Two other weak lines sum up to 6.06 MeV in detectors with moderate energy resolution. A 
yS-electron is emitted on average of 0.3 jis before the 8.78 MeV CK-particle. The copper strip is 
attached to the internal surface of a copper cylinder (not present in Fig. 6 in order to make the 
detector visible) which fully surrounds the detector and serves as a blackbody radiation shield. 
The internal diameter of the cylinder is equal to one of the two copper disks visible in Fig. 6, 
above and below the detector. The copper strip, with an area of about 1 cm^, is placed at the 
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center of the shield. Since the surface area of the main absorber is larger than that of the shields 
in two cases (see Table 2), we expect that a fraction of the CK-particles will hit the absorber 
directly. The yS-electrons, due to their much longer range, can penetrate the shield and deposit 
energy both in the shield and in the main absorber. In the former case, the a events are recognized 
as main-absorber events, while in the latter the /S events are mixed events. 

The detectors were cooled down separately in a low power dilution refrigerator capable of 
reaching a base temperature of -20 mK. The detector thermistors were DC biased through a volt- 
age supply and a load resistor at room temperature between 2 GQ to 20 GD, chosen according to 
the specific sensor characteristics. The typical operating temperature was -25 mK, correspond- 
ing to a thermistor resistance of ~10 MQ. The thermistors used in these tests were optimized for 
temperatures higher than those assumed in the simulations of Section 3. The VRH parameters 
are in this case Rq = 3.1 and Tq = 5.8 K. The static bias voltage across the thermistors ranged 
typically from 10 to 20 mV. Voltage pulses were read out by a DC-coupled low noise diff'eren- 
tial voltage amplifier, followed by a filtering single-ended stage. The front end electronics was 
at room temperature. The signals were acquired by a 12-bit transient recorder, collecting 1024 
points for each pulse, and registered for off'-line analysis. The main purpose of the experiment 
was to verify and understand the surface event discrimination capabilities of the detector and 
energy resolution was not optimized. The data sets discussed in the next section were acquired 
during measurements lasting typically several hours. 

5. Experimental results and data discussion 

The thermal model described in Section 3 for SSBs was used to simulate the detector behavior 
assuming the Cuoricino thermal parameters. In case of energy absorbed in the shield (surface 
event), it predicts a substantial shape and amplitude diff'erence between the signal read out by 
the shield thermistor and that read out by the main-absorber thermistor (see Fig. 3c). A fast and 
high signal in the shield thermistor is accompanied by a simultaneous low and slow signal in the 
main-absorber thermistor. This is the main ingredient which allows the identification of the event 
origin. The first objective of the experimental tests is to verify this behavior, at least qualitatively. 

We have not simulated directly the detectors realized here because ad hoc measurement of 
the thermal parameters have beeen performed only on Cuoricino detectors, and not for the type of 
mounting used in the present set-up. For example, the thermal coupling between main absorber 
and heat sink is diff'erent here from that adopted in Cuoricino, and has never been measured. In 
addition, the heat capacities and operating temperatures for our set-up are very diff'erent from 
those used in Cuoricino. Neverthelss, the basic detector structure is analogous to that simulated, 
and therefore we expect a qualitative agreement between simulations and experimental data. 
This agreement can be clearly appreciated in Fig. 7, which shows experimental pulses from 
both thermistors originated by an a-particle fully absorbed in one of the Ge shields of the SSB. 
The diff'erence between the two pulses is attenuated by a low pass filter used in the experiment 
with a cutoff' frequency immediately below 50 Hz, and limiting therefore the pulse rise time 
to about 6 ms. This filter helps in improving the signal to noise ratio without aff'ecting the 
discrimination capability of the realized SSBs. The pulses acquired with the SSBs provided 
with Ge, Si and Te02 shields exhibit similar features. This is not surprising, since the heat 
capacity of the shield thermistor is of the order of 10"^^ J/K at 25 mK (see Table 1), therefore 
dominating over the heat capacities of the shields independently of the material (see Table 2). 
Moreover, the thermal conductances beteween shields and main absorber are similar, as discussed 
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Figure 7: Experimental pulses generated by an alpha particle fully absorbed in the Ge active shield. 

in Section 4.2. However, this parameter is difficult to reproduce and is probably responsible for 
the observed differences. 

5.7. Event origin identification through pulse-amplitude scatter plots 

In the scatter plot obtained from the Ge-shield detector and shown in Fig. 8a, energy depo- 
sitions in the shield can easily be distinguished from those in the main absorber since they are 
distributed in a band with a steeper slope. The main-absorber event band exhibits a fine structure 
with a sub-band characterized by a slightly higher slope. This eff'ect is due to alpha particles 
impinging on the shield which is not read out. These events produce a thermal response which is 
slightly different from that of those occurring in the main absorber, and therefore they appear as 
a separate population. This eff'ect is explained in more detail for the detector with Te02 shields 
(see Fig. 9 and the related explanation in the text). 

The scatter plot can be compared with the simulation of Fig. 4. The qualitative agreement 
between experiment and simulation is once again evident. The simulation shows a shield pulse 
amplitude saturation which does not appear in the experimental data. The lower temperatures 
used in Cuoricino implies much lower heat capacities with respect to the experimental case. 
This determines a non-linearity of the pulse amplitude-energy relation when particles impact 
the shield, while in the experimental data this relation is linear both for the shield and the main- 
absorber events, since the detectors are operated at 25 mK. This explains why the two main event 
populations in Fig. 8a are distributed along straight lines. 

Bulk events giving pulses in the lower band in Fig. 8a are mainly due to natural y radioac- 
tivity and to cosmic-ray muons passing through the main absorber without interacting with one 
of the shields. The points between the two bands are mainly due to particles which deposit en- 
ergy in both main absorber and shield. Most of these mixed events are due to muons crossing 
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Figure 8: (a) Experimental scatter plot showing the relationship between main-absorber pulse amplitudes and shield 
pulse amplitudes acquired for one of the Ge shields. Two populations of events appear, distributed in bands with different 
slopes. These bands include mainly events from only one of the absorbers, (b) Energy spectrum from the shield read-out: 
the higher slope band in the scatter plot is selected so that only pulses due to total energy depositions in the shield are 
kept. The expected structure of the spectrum of the source, which is outside the detector and faces the shield, is clearly 
extracted from the raw spectrum after surface event discrimination. 



both the read out detector elements, as shown in a previous work [17], and to y's undergoing 
a Compton interaction in one of the two elements followed by a Compton interaction or pho- 
toelectric absorption in the other one. As mentioned in Section 3, the possibility of identifying 
mixed events is important for experiments such as CUORE, as this will add the capability of 
discriminating energy depositions due to Te02 surface contamination, although this contribution 
looks less dangereous than that provided by the inert Cu structure. 

The effectiveness of the discrimination power of this technique is shown in Fig. 8b. The 
structure of the external a- source is clearly extracted from the raw energy spectrum of the shield 
after selecting the higher slope band shown in fig. 8a. 

Data analysis of the Te02 SSB evidenced the possibility to observe energy deposition in one 
of the active shields by reading the thermistor on the other shield. Fig. 9 shows the scatter plot 
obtained from the individual reading of the thermistors on one shield and on the main absorber, 
while acquiring signals not just from these two thermistors, but also from the sensor on the 
second shield. In addition to the usual main-absorber and shield event bands, a second population 
appears in a band with a slope slightly higher than that of the main-absorber band. To identify 
the origin of this population, we selected the shield events band in the scatter plot obtained from 
the thermistor on the other Te02 shield, shown in the inset, and cut on the first graph. The 
result of such selection is that the additional population close to the bulk region identifies pulses 
due to a energy depositions on the other shield. In the shield thermistor read-out, we would, 
in principle, expect no appreciable diff'erence between pulses corresponding to energy releases 
in the other shield or in the main absorber. However, the diff'erent heat flow paths in the two 
cases (an additional thermal impedence is present for energy deposited in the shield) lead to 
slight difl'erences in pulse shapes (see Section 6). Since the pulse amplitudes in the scatter plots 
are determined through a digital optimum filter which uses an average pulse of bulk events as 
signal template, the amplitude reconstruction depends on the pulse shape. This determines the 
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Figure 9: Experimental scatter plot reporting shield pulse amplitudes vs. main-absorber pulse amplitudes for the first 
Te02 shield. In the inset, the same for the second Te02 shield. Surface events due to energy depositions on one shield 
can be reciprocally seen on the scatter plot corresponding to amplitudes seen by the thermometer on the other shield. 



additional population observed in the scatter plot. 

The signals from the shields can be acquired by connecting the thermistors in parallel. This 
interesting possibility would reduce significantly the number of read-out channels, which is a 
critical parameter for cryogenic systems that have to reach temperatures below 10 mK. In fact, 
each read out wire is a source of power load. This possibility was tested with Si and Te02 SSBs, 
and the results from the latter are shown in Fig. 10. The two boxes identify a interactions on 
each shield. The diff'erence in the slope of the bands is due to unequal thermal conductance 
between the two shields and the main absorber, as the same thermal coupling is difficult to 
reproduce experimentally. The use of the Ge stand-offs for the Te02 case (see Fig. 5) has 
however improved considerably the reproducibility of the couplings, making the results on the 
parallel readout particularly significant. 

5.2. Event origin identification through pulse shape analysis 

As mentioned in Section 3, shield events can be distinguished from main-absorber events 
by a pulse- shape analysis of the signal from a shield thermistor, without additional information 
from the main thermistor. This discrimination possibility was confrimed in all the three detectors 
realized, by analyzing in particular the rise time of the pulses from the shield thermistors. 

In order to show the potential of this method, we report the results obtained with the detector 
equipped with Si shields as an example. Fig. 11 shows the rise time vs. pulse amplitude for 
one shield thermistor of the Si SSB. This can be compared to the simulated behavior described in 
the last part of Section 3. The Si shield sees two clearly recognizable groups of pulse amplitudes: 
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Figure 10: Experimental scatter plot for Te02 SSB corresponding to a parallel read-out of the shield thermistors. The 
two boxes contain events observed in the two shields 



fast pulses due to particle interactions in the shield (r^- 10 ms) and slow pulses originating from 
energy depositions in the main absorber (t^~ 30 ms). In the inset, the scatter plot shows that the 
class of fast pulses corresponds to shield events. 

The clear event origin discrimination achievable with shield thermistors offers the possibility 
of reading out the entire detector using a single wire pair. The main-absorber thermistor could 
be eliminated and the shield thermitors could be connected in parallel at the detector level. Only 
one wire pair would be used to extract the common signal from the cryostat. From this signal, 
it would be possible to obtain the fundamental information concerning the energy deposited in 
the main absorber from the analysis of the slow pulses (corresponding to bulk events), while 
the surface events could be cut by rejecting the fast rise time pulses. Also mixed events with a 
substantial energy fraction deposited in the shield would be identified (se Fig. 3(b)). Future tests 
will be needed to measure the eff'ect on the energy resolution relative to the usual main-absorber 
thermistor, and optimize any trade-off' between background rejection and energy resolution acting 
on the thermal parameters of the detector. 

53. Quantification of the surface-event rejection power 

In order to estimate the impact of the proposed method on the sensitivity of future experi- 
ments, it is useful to attempt a preliminary evaluation of the surface-event rejection power. We 
have performed this investigation using the identification method based on the scatter-plot anal- 
ysis. The detector selected for this study is the one with Te02 shields read out in parallel (the 
third one in Fig. 5), since it is the most similar to a practical device for double beta decay search. 
First we have fitted the points corresponding to bulk events with a straight line starting from the 
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Figure 11: Shield pulse rise-time vs. amplitude for the detector with Si shields. Two classes of events corresponding to 
different can be identified. In the inset, scatter plot showing shield-pulse vs. main-absorber pulse amplitudes for one 
Si active shield. The relationship between the two graphs is indicated with arrows: pulses corresponding to shield energy 
depositions can be distinguished because of lower with respect to main-absorber energy depositions. 
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origin in the scatter plot reported in Fig. 10. These points are easily identifyable and correspond 
to the narrow band with the smallest slope and ending at about 2500 mV on the x-axis. Let (3 be 
the slope of this band. We can then construct a parameter, associated to each event and defined 
as SEI (Surface Event Index), in the following way: 

SEl . ^^S^ 1 (5) 

P • Ay {main absorber] 

where (AV [main absorber], AV [shield]) is a pair of values defining a point on the scatter plot 
and corresponding to an event registered by the detector. 

For construction, SEI will be (within experimental errors) for bulk events, while it will be 
> for events characterized by a direct energy deposition in one of the two shields. Looking at 
the distribution of SEI, we expect then a peak centered on and corresponding to bulk events 
and two peaks centered on values higher than and corresponding to shield events. We expect 
also a small population of events lying outside the three main peaks, corresponding to the cases 
of mixed main absorber- shield energy depositions, already discussed in Section 5.1. 

The separation between the bulk event peak and the shield event peaks in the SEI distribution 
allows to evaluate the rejection power of the method. By fitting the peaks with Gaussians, it 
is possible to determine the leak of a Gaussian related to shield events into a region centered 
on and corresponding to bulk events selected with a given efficiency. Of course, we expect 
that the rejection power depends on the threshold for the energy deposited in the main absorber. 
The lower the threshold, the lower the signal-to-noise ratio, with consequent merging of the 
distributions related to bulk and surface events. 

The results of this method applied to the detector with two Te02 shields are appreciable in 
Fig. 12. In the x-axis, we report the pulse amplitude read by the bulk thermistor calibrated in 
energy. The histogram of these pulse amplitudes shows clearly two peaks identifiable as the ^^K 
peak (1460 keV) and the peak (2615 keV) due to the natural y radiactivity in the laboratory. 
This identification makes the energy calibration of the x-axis possible. On the y-axis of the 
main plot of Fig. 12, we report the parameter SEI constructed as described above. As expected, 
it is possible to observe a horizontal band of events at SEI ~ (bulk events), one at SEI ~ 4 
(shield 2 events) and one at at SEI ~ 8 (shield 1 events). In the bands of shield events, the four 
clusters corresponding to the main lines of the a source used to simulate surface contamination 
are appreciable. The inclined shape of these clusters (more evident for events on shield 1) are 
probably explanable as a position eff'ect due to the presence of point-like contacts between the 
shield and the main crystal provided by the Ge stand-off's. If an cK-particle deposits its energy 
close to a stand-off', it generates a higher pulse in the main-absorber thermistor with respect to 
an energy release far from the shield-crystal contact point. The eff'ect is opposite for the shield 
thermistor. 

After selecting an energy window (2000 - 4000 keV) containing the Q- value of ^^^Te Double 
Beta Decay (2527 keV), the three horizontal bands described above are projected into the SEI 
distribution showed in log scale on the left side of Fig. 12. Three peaks appear, corresponding to 
bulk, shield 1 and shield 2 events. The bulk event peak exhibits a tail at high SEI values due to 
the mixed events discussed previously. In order to assess the rejection power, the shield 2 event 
peak has been fitted with a Gaussian, centered on 4.208 and with cr = 0.127. The distribution 
of the bulk events, including the mixed event tail, stops at about SEI = 2.5, about 13 cr's away 
from the centroid of shield 2 events. This analysis shows that in the double beta decay region 
there is no contamination of a's impinging on shield 2 in the bulk event population and that 
the rejection power is essentially 100%. A similar result is achieved when analysing the shield 1 
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Figure 12: The Surface Event Index (SEI) introduced in the text is plotted as a function of the energy deposited in the 
main absorber for the detector with two Te02 shields read out in parallel. The distribution of SEI is projected on the 
left for events in the energy window 2000-4000 keV. Three event populations (bulk, shield 1 and shield 2 events) appear 
clearly. The separation power of the event classes is discussed in the text. 
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events. These preliminary results are very encouraging and indicate clearly that a full recognition 
of surface events is possible with this technology. 

Although the main purpose of the present method is to discriminate surface events emitted by 
Cu surfaces, we believe that it has the potential to reject mixed events as well, especially those 
really dangerous for the search for Ovfip. In Fig. 12, the mixed events are mostly concentrated 
close to the bulk event band. This is due to the fact that in an aboveground experiment most of 
the mixed events are due to energetic cosmic muons with inclined trajectories which cross both 
the main absorber and the shield. Due to the thinness of the shield, these muons normally deposit 
much more energy in the main absorber than in the shield. A fraction of them could anyway be 
discriminated by fitting the SEI distribution around with a symmetric gaussian plus a tail on 
the right side. However, this operation is not relevant here since these events will not be present 
in an underground set-up, where the dangerous background is mainly due to a radioactivity. As 
already pointed out, in order to simulate a Qv^p CK-particles must deposit at least -1.5 MeV in 
the shield. The energy share in this case is such to produce a mixed event similar to that depicted 
in Fig. 3(b), and the corresponding band on the scatter plot is visible in Fig. 4. Therefore, after 
the selection of an energy interval around the OvfiP energy, we expect a SEI distribution with 
a clear gap between the bulk event peak and the mixed event population. In order to quantify 
the rejection power in this case, we are designing an experiment, to be performed underground, 
in which a superficial a source is intentionally deposited on the main-absorber surface, below a 
Te02 shield. 

6. Future prospects: Pulse Shape Discrimination with main-absorber signals 

The prototype SSBs for discriminating surface contamination events described here consist 
of three separate thermistor read-outs, one for the main absorber and one for each of the two 
active shields. Data analysis and simulations show an interesting possibility of doing this dis- 
crimination using only the thermistor connected to the main absorber. This approach obeys to 
a completely diff'erent philosophy with respect to the results presented in the previous Section, 
since the concept of auxiliary bolometers is abandoned and the shields play a passive role aiming 
at the modification of the signal shape. This method looks very promising for real application in 
future experiments, and will be discussed with some detail in this Section. 

Fig. 13 shows the distribution of decay time (90%-30%) vs. amplitude for pulses from the 
main absorber of the Te02 SSB. The plot shows that the typical decay time is of the order of 
22 ms, rather independent of amplitude above a given threshold. The spread of the decay times 
at small amplitudes is an eff'ect ascribable to noise and especially to 'pile-up' pulses. These latter 
are a class of partially-overlapping pulses (in the time domain) whose tricky superposition leads 
them to be recognized as a single pulse with higher decay time. The aboveground operation of 
the detector and the intrinsic slowness of bolometers make the pile-up eff'ect particularly severe. 
The plot shows also diff'erent groups of pulse amplitudes corresponding to longer decay times. 
These groups are included in the box outlined by dotted lines and correspond to the points circled 
in the inset, which shows the familiar scatter plot of shield pulse amplitudes vs. main-absorber 
pulse amplitudes. This shows clearly that shield events are identifiable as longer decay time 
pulses when read out by the thermistor on the main absorber. The additional population selected 
close to the bulk event band is due to events occuring in the other shield and discussed previously 
(see Fig. 9 and related discussion). The event discrimination in Fig. 13 using decay times is not 
as clear as in the inset, however further optimization of the thermal parameters may be possible 
in the future. 
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Figure 13: Main-absorber pulse decay time vs. pulse amplitude for the detector with Te02 shields. A structure with 
groups of events with higher rj can be isolated from the region that identifies usual bulk events, characterized by ~ 22 
ms decay time rather independent of the amplitude. In the inset, scatter plot of shield-pulse vs. main-absorber pulse 
amplitudes using one Te02 shield. As can be seen in the two graphs, events in the shields can be separated from the main 
absorber events by selecting events with high rj. 
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In order to understand better this effect, the decay time as a function of the pulse ampHtude 
was analyzed through our pulse shape simulation tool, based on the thermal model of Section 
3. The simulations confirm that a longer decay time is expected for events originating in the 
shield. This is due to the additional thermal impedance between shield and main absorber, which 
slows down the thermal pulses as seen by the main-absorber thermistor. An important point 
emerging from the simulation results is that this thermal eff'ect is enhanced by the presence of 
an element with large heat capacity connected to the shield. In ordinary conditions, this element 
is the electron system of the shield thermistor. Another result is that the eff'ect does not require 
the shield thermistor to be linked to the heat bath through the wires for signal acquisition. In 
principle, therefore, the thermometer can be replaced by any passive element with appropriate 
heat capacity with the purpose to enhance the pulse shape diff'erence. 

To analyze quantitatively this possibility, a series of simulations were carried out with a 
5x5x5 cm^ Te02 cubic main absorber and a single Si shield with a small Cu block thermally 
attached to the Si. In this model there is only one thermistor attached to the main absorber, and 
none on the shield. The simulations estimate the decay time vs. amplitude for an energy range 
100 keV - 9 MeV. The size of the Cu block was varied, spanning a heat capacity range from 
6.1xlO"^-T J/K to 4.9xlO"^-T J/K, where T is measured in Kelvin. The base temperature was 
set at 9 mK, and all other parameters listed in Table 1 were used, in order to simulate Cuoricino 
/ CUORE type detectors. The results of these simulations for several diff'erent heat capacities of 
the Cu block are illustrated in Fig. 14. It shows the ratio of the decay-time of a shield event 
to a main-absorber event as a function of the pulse amplitude, as observed by the main-absorber 
thermistor. The simulations show clearly the increasing trend of this ratio as the size of the Cu 
block increases. This technique is potentially a very powerful tool for large-scale experiments 
with a large number of detectors. It allows for shield event discrimination resulting from surface 
contamination thanks to pulse shape analysis techniques and without the complexity and heat 
load of additional thermistors and their accompanying read-out wires. 

A price to pay is a reduction of the pulse amplitude, as clear from Fig. 14 when one looks 
at the maximum amplitude corresponding to 9 MeV energy deposition. However, a reasonable 
compromise can be found. For example, the second smallest copper block (which has a heat 
capacity at -10 mK -200 times higher than that of the Si shield and similar to that of the ther- 
mistor) determines an acceptable pulse amplitude reduction of the order of 10%, providing a 
decay time longer by the same amount. With the very high signal-to-noise ratio pulses expected 
in the OvyS/? region, this diff'erence in shape should be easily identified. The capability fo this 
method to identify also mixed events needs to be investigated. We expect that if a considerable 
amount of the energy is deposited in the shield the slow time constant related to the discharge 
of this energy into the main absorber should modify the pulse shape read by the main-absorber 
thermistor. An experimental confirmation is however mandatory. 

In order to simplify the detector structure and to avoid the Cu block, one could choose a shield 
material with a high specific heat at low temperatures. This material cannot be a normal metal, 
for which the specific heat is proportional to the temperature and definitely too high, given the 
large size of the shield. One should therefore choose a dielectric or a superconductive material 
(in which the electron contribution to the specific heat vanishes exponentially). However, this 
choice is not straight-forward, as this material should satisfy many constraints simultaneously: to 
have a low Debye temperature in order to get sufficiently high heat capacities; to be easily put in 
the form of large surface single crystals; to exhibit a very low intrinsic bulk radioactivity; to have 
surfaces easily cleanable by polishing and / or etching; to have thermal contractions which match 
reasonably well the Te02 ones. This means to open a new R&D line with uncertain outcome. 
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Figure 14: Ratio between the decay time rj for main-absorber pulses corresponding to shield energy deposition and 
that corresponding to main-absorber energy depositions vs. main-absorber pulse amplitudes. The range of the energy 
deposited in the main absorber is 0-9 MeV. The simulated detector is a S SB with a Te02 cubic main absorber (5 cm 
side length), a single Si active shield (5x5 cm^ surface and 300 fim thickness) and a small Cu block with variable size 
(maximum volume 0.5 mm^) glued on the shield. The graph shows that the Td ratio increases with the heat capacity of 
the Cu block. 
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Seeking simplification, another solution could be to decrease substantially the thermal conduc- 
tance between shield and main absorber, in order to enhance the pulse shape diff'erence without 
introducing an additional heat capacity. This is not easy too, since the mechanical fixing of the 
shield requires three - four contact points with surface of the order of a square millimeter (a 
smaller area would not provide a safe attachment). Since at low temperatures the thermal impe- 
dences are dominated by the contacts between dissimilar materials rather than by the conduction 
along the bulk, the order of magnitude of the thermal conductance is fixed by the geometry. These 
considerations show that the Cu block solution looks the most viable, providing in addition an 
easy method to tune the pulse modification eff'ects just acting on the block volume. 

7. Conclusions and future prospects 

In conclusion, the ability to discriminate surface events from bulk events has been demon- 
strated experimentally and a thermal model describing the system has been developed. The 
model agrees with experimental measurements qualitatively and future application for bolometric 
experiments seems promising. Shield events can be distinguished from absorber events through 
amplitude comparison of the simultaneous pulses in the main absorber and the active shields. 
Reading out the shield thermistors in parallel is an eff'ective way of reducing the total number 
of acquisition channels without compromising the event location discrimination of the SSBs, 
provided that reasonable reproducibility in the thermal couplings is achieved. 

Pulse shape discrimination can help discriminate main absorber from background events and 
reduce the overall number of acquisition channels further: two diff'erent ways could be explored. 
Event location can be identified through pulse shape discrimination of the active shields read 
out in parallel, in principle eliminating the needs for the main-absorber read-out. This method 
can diff'erentiate between shield and mixed events from main absorber ones through the diff'er- 
ent resultant rise times, and would provide information on the energy deposited in the main 
absorber through the slow pulses; it probably reduces the overall energy resolution of main ab- 
sorber events, but is worth being studied and optimized. The other way, acquiring only the signal 
coming from the main thermistor, seems promising as well by diff'erentiating by the event de- 
cay time. The shields would act as pulse- shape modifiers, without the complication of additional 
read-out channels from independent thermistors. The addition of a thermal load on the shield can 
enhance this eff'ect. This will be very useful in large-scale experiments for which surface-event 
discrimination is necessary as much as the lowering of the overall number of readout channels. 

Active discrimination of surface background may have a large impact on the sensitivity of 
next-generation CUORE-like double beta decay experiments. For the CUORE experiment, the 
background in the OvySyS region has been extensively studied with Monte Carlo simulations, taking 
into account all anticipated sources (bulk and surface radioactive impurities, cosmic-ray muons, 
neutrons, etc.) [7]. The results show that the present configuration with no SSB has a background 
level of about 10 counts/keV/ton/y, providing a sensitivity to the neutrinoless double-beta-decay 
half-life of ^^^Te of about 2x10^^ y. This sensitivity would just reach the inverted hierarchy re- 
gion of the neutrino mass pattern [30, 31, 32]. By eliminating surface events, an additional factor 
of 10 reduction in background can be achieved according to simulations. In terms of sensitivity 
to Majorana neutrino mass, this would extend the search deep into the inverted hierarchy region 
and substantially increase the discovery potential of the experiment. 
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